Computational study of the dielectric properties of [La,Sc]203 solid solutions 
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First-principles calculations were used to compute the dielectric permittivities of hypothetical 
[La,Sc]203 solid solutions in the cubic (bixbyite) and hexagonal La203 phases. Dielectric enhance- 
ment is predicted at small Sc concentrations due to the rattling ion effect. Similar calculations for 
a model amorphous La2 03 structure show little change in permittivity when a small amount of 
Sc is substituted for La. In this case, the local environment around the Sc changes in a way that 
compensates for the rattling ion effect. 



High permittivity ("high k") dielectric materials 
for microelectronics applications require a combina- 
tion of properties, including a permittivity k larger 
than that of amorphous Si02 (k = 3.9), thermody- 
namic stability in contact with Si, and sufficient band 
offsets with respect to Siji While high-K dielectrics 
based on Hf02 have been incorporated into commer- 
cial computer chips, ^ other "higher-K" materials con- 
tinue to be actively investigatedf^ as future Metal-Oxide- 
Semiconductor Field-Effect Transistor (MOSFET) scal- 
ing will likely require dielectrics with higher k than 
Hf02.^ An additional concern for dielectrics for MOS- 
FET applications is whether the dielectric is amorphous. 
Amorphous dielectrics are considered superior, as micro- 
crystalline dielectrics can exhibit significant charge leak- 
age and ionic diffusion due to grain boundaries.^ On the 
other hand, crystalline dielectrics are still of interest, as 
epitaxial dielectric growth^ would avoid the above prob- 
lems. 

Early compilations of candidate high-K materials in- 
dicated that La203 has higher permittivity than Hf02,^ 
thus much attention on higher-K materials has focused on 
La203 and other i?203-based materials (i? = rare earth). 
High permittivities have been reported in both single 
component i?203 compounds, such as Sm203,^ and in 
stoichiometric i?Sc03 compounds, such as GdSc03 and 
DySc03.-' This work explores a possible way to fur- 
ther enhance the permittivity of rare-earth dielectrics: 
nonstoichiometric substitution of Sc for the R ion. 
This idea is inspired by recent theoretical work on the 
Bai_a;Caj;Zr03 (BCZO) perovskite system.'^ BCZO has 
a dielectric anomaly, where the permittivity is higher for 
intermediate compositions than for either endmember, 
reaching a maximum in the range x « 0.1 to 0.2i^i^ The 
enhancement in permittivity is due to a "rattling ion" ef- 
fect: the smaller Ca ion in the lattice site of the larger Ba 
ion is relatively loosely bound and can thus move farther 
under an applied electric field. A similar effect might be 
expected when a small Sc ion substitutes isovalently for a 
larger R ion in R2O3 compounds. The largest R^~^ ion is 
La'^+ (Ref. [13) ; thus we explore the dielectric properties 



of the La2_22;Sc22;03 (LSO) solid solution system in this 
work. 

La203 has two room-temperature polymorphs: a cubic 
(c) bixbyite phase and a hexagonal (h) phasefi2,iii Due 
to the small difference in energy between c-La203 and 
/i-La203 (which is thermodynamically stable at sinter- 
ing temperatures of order 1000 °C), it is not clear which 
phase is the most stable at room temperature.^'' Wc thus 
explore both structures. SC2O3 has the cubic bixbyite 
structure at room temperature. LSO has one ordered in- 
termediate phase, with the perovskite {p) structure>i^ 

We ran first principles calculations on the 40-atom 
primitive cell for c-La203 cell, and on a 40 atom /i-La203 
supercell (each cell parameter doubled), for both pure 
La203 and for the same cells with one Sc substituted for 
La {x — 0.0625). The initial atomic positions were taken 
from the literatureiii Additional calculations were per- 
formed for: (1) c-LSO with x — 0.25, to test the linearity 
of K vs. x; (2) c-LSO with x ^ 0.9375 ("reverse" Sc 
— La substitution); (3) ideal p-LaSc03; (4) a model for 
amorphous La203, and (5) a model for amorphous LSO 
(a-LSO) with x = 0.0625. All calculations were per- 
formed at the experimental volume or the volume given 
by Vegard's law: V = {1 - x)Vl!,^o^ -\- x VscaOa- 

Structural relaxations and dielectric constant calcula- 




FIG. 1: (Color online) Relaxed models for (a) amorphous 
La203 and (b) La2-2:.Sc2:r03 {x = 0.0625). 
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FIG. 2: (Color online) Comparative computed dielectric con- 
stant for La2-2a;Sc22!03 Structures studied in this work. 
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tions were perforraed using the density functional theory 
code PHASE,— Ultrasoft pseudopotentials^ were used, 
with the PBE generalized gradient approximation^^ for 
the exchange-correlation functional, and a plane wave 
basis set with a 36 Rydberg cutoff energy for the elec- 
tronic wavefunctions. Brillouin zone integrations were 
performed with the F point only for the force calcula- 
tions, and were done with appropriate grids for the elec- 
tronic dielectric constant and the Berry phase calcula- 
tions. All structures were relaxed until the force on each 
atom was smaller than 0.002 Ry/Bohr. For c-La203, 
there are two inequivalent La positions: 8(6) and 24(d) 
(Ref. [Ill; modern convention for space group /a3 Wyck- 
off positions from ReL 16^. At both x = 0.0625 and 
X = 0.9375, the lowest energy is found when the minor- 
ity species is on the 8(6) site. Because a single occupied 
8(6) is lower in energy than 24((i), it seems reasonable 
to assume that, for lowest energy configurations, the 8(6) 
sites fill with Sc before the 24(d) sites. For the x — 0.25 
calculation, therefore it was assumed that lowest energy 
occurs when the 8(6) site is fully occupied. 

A model for the La2 03 amorphous phase was gen- 
erated by the ab initio molecular dynamics "melt and 
quench" techniquCflS A 80-atom cubic La203 cell was 
"melted" at 3000 K until the root mean square displace- 
ment of oxygens from their original positions was greater 
than the original oxygen-oxygen near-neighbor distance. 
Then the cell was quenched for 1 ps of simulation time at 
1000 K. The lowest energy structure during this run was 
recorded, and subsequently a full relaxation (cell shape 
and atomic positions) was performed to obtain an a- 
La2 03 model. (The melt-and-quench technique was done 
using the VASP package^ with projected augmented 
wave functions^^ and the local density approximation for 
exchange-correlation; the final relaxation and subsequent 
steps used PHASE, as above.) After this relaxation, 
bond valence sums^i^ were computed to identify the 
most underbonded La as the most favorable positions for 
Sc. The Sc was placed on the most underbonded La site 
because the smaller ionic radius of Sc favors the lower co- 



ordination number. The two most underbonded La sep- 
arated by more than a/2 were replaced by Sc, and the 
structure then fully relaxed under the new cell volume to 
obtain a model for a-LSO; x = 0.0625. The procedure 
insures that the a-LSO model is as similar as possible to 
a-La203 model for direct comparison of their properties 
(Fig. HI). ^ 

The dielectric properties were computed using the 
frozen phonon method, as discussed in Ref. ^2^. Briefly, 
the static dielectric constant of a crystal is given by 

K = K<'l + Ki°" = K°' + C^^, (1) 

where C is a constant, k*'' the electronic contribution to 
the dielectric constant, k'™ the ionic contribution, 
the dipole moment associated with mode /i and lu^ its 
frequency. Aside from obtaining the dielectric constants, 
the breakdown of the lattice dielectric contribution over 
individual modes allows one to find the origins of any 
dielectric anomaly that is due to lattice effects. A rattling 
ion effect will show up as a significant contribution to k 
from a mode with low frequency and/or high that is 
located on or near the rattling ion. 

Fig. [5] shows the dielectric constant computed for each 
structure. Table |T] shows, in addition, k'™, and the 
calculated bandgaps. For small x, the dielectric constants 
of both c-LSO and h-LSO are higher than the weighted 
average of the endmember La203 and SC2O3 compounds. 
Phonon density of states calculations (Fig.|31) show that, 
as x increases, the frequencies of the modes around 200 
cm-i to 300 cm that contribute significantly to the 
dielectric constant decrease. The dielectric constant is 
thus enhanced, per Eq. ([T])i^ Analysis of these lower- 
frequency phonons shows that the eigenvectors are dom- 
inated by Sc motion coupled to motion of its nearest 
neighbor O octahedron in the opposite direction. We 
therefore conclude that the enhancement of k is due to 
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FIG. 3: (a) Comparative phonon density of states for c- 
La2-2a;Sc2i03 versus X. (b) Contribution of phonons to di- 
electric constant, versus x, binned by 10 cm~^ range of fre- 
quency. 



the rattling ion effect. This hypothetical effect is distinct 
from the enhancement of permittivity due to rare-earth 
doping of Hf02 , which is believed to result from the sta- 
bilization of the Hf02 tetragonal phase relative to the 
monoclinic one.-^- For the /i-LSO case, the substituted 
Sc atom displaces along the hexagonal c axis to optimize 
the structure. Tight (0.201 nm) and weak (0.304 nm) 
Sc-0 bonds are formed. As a result, Sc moves more eas- 
ily along the c axis, but is more rigid in the ab plane. 
The Kcc component is enhanced, while Kaa and Kbb are 
reduced. Overall, k is enhanced. When a small amount 
of La is substituted for Sc in C-SC2O3, the k^^ increases, 
which we attribute to the bandgap reduction (Table |I]) , 
but was not calculated, as a mode with imaginary 
frequency {i.e. an instability) appears. In other words, a 
small amount of La substitution is predicted to break the 
cubic symmetry. The relative increase in k with compo- 
sition is less than that found for BCZOj>^ and the dielec- 
tric constant predicted at 2: = 0.25 does not exceed that 



computed for the intermediate perovskite LaScOa phase. 
Nonetheless, if a phase that is close in composition to 
i?203 is desired, our results suggest a way to increase the 
permittivity with relatively small doping of Sc. Further- 
more, if it were possible to substitute a small amount of 
Sc on the A site of a i?Sc03 perovskite (an isostructural 
analog to BCZO) without destabilizing the perovskite 
phase, the rattling ion effect might yield a permittivity 
even higher than for stoichiometric p-RScOs, compounds. 

Our results predict that k for amorphous La203 is 
higher than for c-La203 but smaller than for ft,-La203. 
Delugas et al. predict a similar permittivity enhance- 
ment for amorphous SC2O3 with respect to c-Sc203i^ In 
contrast to the crystalline cases studied, the computed 
dielectric constant of a-La203 does not increase when a 
small amount of Sc is substituted for La. We explain the 
lack of any rattling ion effect as a consequence of different 
Sc and La local environments in a-LSO. Upon relaxation, 
the Sc ions have on average 5.0 oxygen neighbors at an 
average distance of 0.208 nm, while the La ions have on 
average 6.0 oxygen neighbors at an average distance of 
0.246 nm. The results suggest that the local environment 
around Sc in a-LSO "tightens up" relative to the typical 
La local environment in a way that compensates for the 
rattling ion effect. 

To summarize, small substitution of Sc for La in crys- 
talline La203 compounds is predicted to cause an anoma- 
lous increase in the dielectric constant due to the rattling 
ion effect: a small Sc ion in a loose environment can move 
farther in an applied electric field. For amorphous La203, 
we predict no anomalous increase, as La and Sc have dif- 
ferent local environments in this case. 
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